Mutually unbiased bases, orthogonal Latin squares, and hidden-variable models 
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Mutually unbiased bases encapsulate the concept of complementarity — the impossibility of si- 
multaneous knowledge of certain observables — in the formalism of quantum theory. Although 
this concept is at the heart of quantum mechanics, the number of these bases is unknown except 
for systems of dimension being a power of a prime. We develop the relation between this physical 
problem and the mathematical problem of finding the number of mutually orthogonal Latin squares. 
We derive in a simple way all known results about the unbiased bases, find their lower number, and 
disprove the existence of certain forms of the bases in dimensions different than power of a prime. 
Using the Latin squares, we construct hidden-variable models which efficiently simulate results of 
complementary quantum measurements. 

PACS numbers: 03.65.Ta, 02.10.Ox 
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I. INTRODUCTION 

Complementarity is a fundamental principle of quan- 
tum physics which forbids simultaneous knowledge of cer- 
tain observables. It is manifested already for the simplest 
quantum mechanical system — spin- ^ ■ If the system is 
in a definite state of, say, spin along x, the spin along y 
or z is completely unknown, i.e., the outcomes "spin up" 
and "spin down" occur with the same probability. The 
eigenbases of a x , <J y and <j z Pauli operators form so-called 
mutually unbiased bases (MUBs) : Every vector from one 
basis has equal overlap with all the vectors from other 
bases. MUBs encapsulate the concept of complementar- 
ity in the quantum formalism. Although complementar- 
ity is at the heart of quantum physics, the question about 
the number of MUBs remains unanswered. Apart from 
being of foundational interest, MUBs find applications in 

Suantum state tomography [l[ , quantum-key distribution 
| , and the mean King problem [3J . 
A d-level quantum system can have at most d + 1 
MUBs, and such a set is referred to as the complete set 
of MUBs. In 1981 Ivanovic proved by construction that 
there are indeed d + 1 complementary measurements for 
d being a prime number [J]. This result was general- 
ized by Wootters and Fields to cover powers of primes 
[l(. For other dimensions the number of MUBs is un- 
known, the simplest case being dimension six. A consid- 
erable amount of work was done towards understanding 
this problem. New proofs of previous results were estab- 
lished p, [y, 0, Q and t he p roblem was linked with other 
unsolved problems [9|, LL(|- It was also noticed that it 
is similar in spirit to certain problems in combinatorics 
[ll|, [Ej LL3| and finite geometry [1J|. Here, we build upon 
these relations. 

We describe the problem of the number of orthogonal 
Latin squares (OLSs), which was initiated by Euler [15| 
and still attracts lots of attention in mathematics. Al- 
though this problem is not solved yet in full generality, 
more is known about it than about the number of MUBs. 



Using a black box which physically encodes information 
contained in a Latin square, we link every OLS of order 
being a power of a prime with a MUB. For dimension 
six, our method gives three MUBs, which is the maximal 
number found by the numerical research [lfj, [ll| . Utiliz- 
ing known results for OLSs we derive a minimal number 
of MUBs, and disprove the existence of certain forms of 
MUBs for arbitrary d. Finally, using OLSs we construct 
hidden-variable models that efficiently simulate comple- 
mentary quantum measurements. 



II. ORTHOGONAL LATIN SQUARES 

A Latin square of order d is an array of numbers 
{0, ...,d— 1} where every row and every column con- 
tains each number exactly once. Two Latin squares, 
A = [Aij] and B = [By], are orthogonal if all ordered 
pairs (Aij,Bij) are distinct. There are at most d — 1 
OLSs and this set is called complete. The existence of L 
OLSs is equivalent to the existence of a combinatorial de- 
sign called a net with L + 2 rows [16j • The net design has 
a form of a table in which every row contains d 2 distinct 
numbers. They are split into d cells of d numbers each, 
in such a way that the numbers of any cell in a given row 
are distributed among all cells of any other row. The ad- 
ditional two rows of the net correspond to orthogonal but 
not Latin squares, with the entries Aij = j and Aij = i. 

The following algorithm allows us to construct the net 
from a set of OLSs: 

(i) Write the squares in the standard form in which the 
numbers of the first column are in ascending order (by 
permuting the entries, it is always possible to write the 
set of OLSs in the standard form without compromising 
Latiness and orthogonality). 

(ii) Augment the set of OLSs by the two orthogonal 
non-Latin squares j4y = j and Ay = i. 

(iii) Write the rows of the squares as cells in a single 
row of the table. The number of the table's rows is now 



equal to the number of squares in the augmented set. 

(iv) In the row of the table which corresponds to the 
square Aij — j, referred to as the "coordinate row," re- 
place the number A^ in the ith cell by A' i4 = id + j, 
where d is the order of the square. 

(v) In every cell of the other rows replace number Bij 
on position j by the integer associated to the number Bij 



B, 



BL 



of the jth cell in the coordinate row, i.e. 
jd + Bij. 

We shall prove that the table generated by this pro- 
cedure is indeed a net design. We use another property 
defining the design: Two numbers in one cell do not re- 
peat in any other cell. This already includes that any 
two cells of two different rows share exactly one common 
number, as if there were no common numbers shared by 
these cells, there would have to be at least two common 
numbers shared by other cells. 

Due to the definitions of A' i4 and B'- 4 and the fact that 
the columns of B^ contain all distinct numbers 0, ..., d—1, 
every row of the table contains d 2 distinct numbers 
0, ..., d 2 — 1. By construction, the numbers of any cell of 
the coordinate row are distributed among all the cells of 
all the other rows. Therefore, it is sufficient to prove the 
property of the net for the remaining rows. Assume to the 
contrary, that two numbers repeat in two cells of different 
rows, (jd + B^ , j'd + Bij> ) = (ld+ Cki , I'd + Cw)- Since 
€ {0, ...,d — 1} the equality can only 
C'kj and Biji — Ckj', i.e., there are rows 
of the squares B and C which contain the same num- 
bers, in the columns defined by j and j' . This, however, 
cannot be because one can always permute the entries 
of, say, square C such that its fcth row becomes the ith 
row (without compromising orthogonality) and the two 
squares would not be orthogonal. 



III. QUBIT 

Consider the squares for d — 2. We link them with 
the complementary measurements of a qubit. The aug- 
mented set of orthogonal squares reads as 



3,3 



, t, i , -Dij j ^ij 



hold if Bij = 
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1 1 



1 

1 0. 



(1) 



The right-hand side square is Latin, the left and middle 
square are orthogonal to each other and to the Latin 
square. These three squares lead to the following net 
design on the left-hand side, in which the numbers are 
represented by pairs mn in modulo-two decomposition: 



6 = 


6=1 


00 01 


10 11 


00 10 


01 11 


00 11 


01 10 



6? 



6? 



6? 



(2) 



On the right-hand side, we write down the complemen- 
tary questions associated with each row. They are an- 
swered by pairs mn in the left- and right-hand column 



of the net design (left column — * answer 0, right column 
— > answer 1). In this way, the questions are linked to the 
orthogonal squares. 

The complementary questions can be answered in 
quantum experiments involving MUBs. Consider a de- 
vice encoding parameters to and n via application of the 
unitary U — cr™<7™. When it acts on |z±) states, they get 
a phase dependent on n and are flipped m times. Thus, 
knowing the initial state, a final measurement in the <x z 
eigenbasis reveals to, giving the answer to the first com- 
plementary question. Similarly, taking |x±) and |j/±) as 
initial states, the results of a x and a y measurement an- 
swer the second and the third complementary question, 
respectively. 



IV. PRIME DIMENSIONS 

For prime d the net has d + 1 rows. The entries of the 
rows corresponding to the OLSs are generated from the 
following formula: 

n = am + 6, (3) 

where the integer a = 1, ..., d — 1 enumerates the rows of 
the table, while the integer 6 = 0, ...,d— 1 enumerates 
different columns, and the sum is modulo d. Additional 
two rows correspond to the questions about to and n, 
respectively. The table for the rows corresponding to the 
OLSs is built in the following way: 

(i) Choose a row, a, and the column, 6. 

(ii) Vary m = 0, ..., d — 1 and compute n using ((3]). 

(iii) Write pairs mn in the cell. 

For example, for d = 3, one has 



6 = 


6= 1 


6 = 2 


00 01 02 


10 11 12 


20 21 22 


00 10 20 


01 11 21 


02 12 22 


00 11 22 


01 12 20 


02 10 21 


00 12 21 


01 10 22 


02 11 20 



TO = 6? 


n = 6? 


n 


= TO + 6? 


n 


= 2to + 6? 



(4) 



The complementary questions are given on the right- 
hand side. Different values of 6 enumerate possible an- 
swers. 

We shall see, again, that the complementary questions 
can be answered using MUBs. Consider encoding of pa- 
rameters to and n via application of U = X m Z n , where 
the Weyl-Schwinger operators X m Z n span a unitary op- 
erator basis. In the basis of Z, denoted as \n), the two 
elementary operators satisfy 



Z\K)=Tft\K) 



X\k) = \k+1), 



(5) 



where r)d = cxp (i2ir/d) is a complex dth root of unity. 
For the same reasons as for a qubit, the first two questions 
are answered by applying U on an eigenstates of Z and 
X operators, and then by measuring the emerging state 
in these bases. 



In all other cases the action of the device is U = 
X m Z am+b = X m Z am Z b . The elementary operators do 
not commute, instead one has ZX = rjdXZ, and it fol- 



the following net design: 



lows that X m Z a 



''/,/ 



- -kamim 



(m-l) 



(XZ a ) m . Finally, the 



action of the device is, up to the global phase, given 
by U ex (XZ a ) m Z b . The eigenstates of the XZ a op- 
erator, expressed in the Z basis, are given by \j) a — 

(i/Vd)J2i= l Vd : ' K ' aSK \ K )^ whcrc s « =„ K + •■• + ( d ~ l ) 

@, and the Z operator shifts them, Z\j) a = \j — 1) . 
After the device, \j) a is shifted exactly b times and sub- 
sequent measurement in this basis reveals the answer to 
the ath question. On the other hand, the eigenbases of 
XZ a for a = 1, ..., d — 1 and eigenbases of X and Z are 
known to form a complete set of MUBs Q . Not only the 
number of MUBs is the same as the number of OLSs, but 
they are indexed by the same variable, a. This allows to 
associate MUB to every OLS for prime d. 



V. POWERS OF PRIMES 

If d is a power of a prime, a complete set of OLSs is 
obtained using operations in the finite field of d elements, 
and one expects that a complete set of MUBs also follows 
from the existence of the field. Indeed, explicit formulae 
for MUBs in terms of the field operations were presented 
in |l|, |7|, |8|. Here, we prove this result in a simple way 
related to [17|, using the theorem of Bandyopadhyay et 
al. d, [lj| : If there is a set of orthogonal unitary matrices, 
which can be partitioned into M subsets of d commuting 
operators, then there are at least M MUBs. They are 
the joint eigenbases of the d commuting operators. 

To illustrate the idea, consider again prime d. Take the 
orthogonal unitary operators S mn = X rn Z n with their 
powers mn taken from the first column of the net. The 
cell of the first and second row corresponds to the eigen- 
bases of Z and X, respectively, whereas the other two 
rows are defined by b — 0, i.e., n = am. According to the 
commutation rule of the elementary operators X and Z, 
S mn and S m ' n i commute if and only if mn' — m'n = 
mod d. Thus, for a fixed row, i.e., fixed a, the set of d 
operators S mn commute, because m{am') — m'(am) = 0, 
and, due to the mentioned theorem, there is a set of d+ 1 
MUBs. 

For d = p r being a power of a prime, the OLSs and 
the net are generated by the formula 



n = a mi 



(6) 



where and © denote multiplication and addition in 
the field, a, 6, m,n £ F^ are field elements, and a =/= 0. 
The first two rows of the table are defined by m = b 
and n — b. In the case of d = 4, the four elements 
10, 1, u), u) + 1} of the field F 4 {ui is the root of x 2 + x + 1 
|17|). when indexed with the numbers {0, 1, 2, 3}, lead to 



00 01 02 03 


10 11 12 13 


20 21 22 23 


30 31 32 33 


00 10 20 30 


01 11 21 31 


02 12 22 32 


03 13 23 33 


00 11 22 33 


01 10 23 32 


02 13 20 31 


03 12 21 30 


00 12 23 31 


01 13 22 30 


02 10 21 33 


03 11 20 32 


00 13 21 32 


01 12 20 33 


02 11 23 30 


03 10 22 31 



(7) 

We use the concept of a basis in the finite field F^. It 
consists of r elements ej, with i = 1, ...,r. Every basis has 
a unique dual basis, ej, such that tr(ej0ej) = Sij, where 
the trace in the field, tr(x), maps elements of F^ into the 
elements of the prime field ¥ p . We use lowercase tr for the 
trace in the field in order to distinguish it from the usual 
trace over an operator, which we denote by Tr. It has 
the following useful properties: tr(ir (By) — tr(x) + tr(y), 
and tr(a x) = atr(x), where operations on the right- 
hand side are modulo p and a is in the prime field. We 
decompose m in the basis e,-, m = m\ e\ © ... © m r e r , 
where rrii = tr(m ej), and n in the dual basis, n — 
n\ ei © ... © n r e r , with m = tr(n ej). Due to the 
properties of the trace in the field and the dual basis 



tr(m n) — > viiirn — fa ■ ft, 



=i 



(8) 



where m — (mi, ...,m r ) and ft — (m, ...,n r ) have com- 
ponents in the prime field, i.e., numbers {0, ...,p— 1}. 



Consider operators defined by the decomposition of m 
and n, S, m = XpZ^ ... X^Z^, where, e.g., Xp 
is the unitary operator acting on the ith p-dimensional 
subspace of the global <i-dimensional space. Operators 
Stfin form an orthogonal basis. They commute, if and 
only if m • ft' — fh! ■ n = mod p. Take the operators 
corresponding to a fixed row of the first column of the 
net, i.e., a is fixed, b = and therefore n — a m. From 
Eq. (|SJ), all these d operators commute if tr(m0o0m ; ) = 
tr(ro' 0a0m), which is satisfied due to associativity and 
commutativity of multiplication in the field. Therefore, 
their eigenbases define MUBs. Again, each row of the 
table is linked with the MUB. 

To make an illustration, consider again the example 
of d = 4. Choose (ei,ea) = (u>, 1) as a basis in the 
field, such that the numbers m are decomposed into pairs 
m — > mi m-2 in the usual way: 0— >0 0, 1 — ► 1, 2 — ► 1 0, 
3 — ► 1 1. The dual basis reads as (e"i,e2) = (l,w + 1), 
which implies that the numbers n are decomposed into 
pairs n — ► m n 2 as follows: 0— >00, 1— >10, 2— » 11, 
3 — > 1. Each pair of numbers of table ([7]) is now written 



vertically as a combination of two pairs of numbers: 



00 01 01 00 


00 01 01 00 


10 11 11 10 


10 11 11 10 


00 00 01 01 


10 10 11 11 


00 00 01 01 


10 10 11 11 


00 00 10 10 


01 01 11 11 


01 01 11 11 


00 00 10 10 


00 10 00 10 


00 10 00 10 


01 11 01 11 


01 11 01 11 


00 01 11 10 


01 00 10 11 


01 00 10 11 


00 01 11 10 


00 10 01 11 


00 10 01 11 


01 11 00 10 


01 11 00 10 


00 01 10 11 


01 00 11 10 


01 00 11 10 


00 01 10 11 


00 11 01 10 


00 11 01 10 


01 10 00 11 


01 10 00 11 


00 00 11 11 


01 01 10 10 


01 01 10 10 


00 00 11 11 


00 11 00 11 


00 11 00 11 


01 10 01 10 


01 10 01 10 



(9) 
MUBs are formed by the eigenbases of operators 
a™ 1 a™ 1 £g> <r™ 2 <r™ 2 , where the powers are taken from the 
first column of this table. The result is in agreement 
with other methods [a, [6[ . The complementary questions 
answered by the states of these MUBs are formulated 
in terms of individual bits mi, TO2, n\, 772, which are 
encoded by U = a™ 1 a™ 1 g) cr™ 2 <x™ 2 . For example, the 
question of the last row is about the values of mi + ri\ 
and 777,2 + "2 ■ 

An interesting feature strengthening the link between 
MUBs and OLSs is the existence of the set of OLSs and 
MUBs which cannot be completed. For example, the 
following net design 



00 01 02 03 


10 11 12 13 


20 21 22 23 


30 31 32 33 


00 10 20 30 


01 11 21 31 


02 12 22 32 


03 13 23 33 


00 11 22 33 


01 12 23 30 


02 13 20 31 


03 10 21 32 



(10) 
cannot have more rows. The MUBs related to this table 
are the eigenbases of X, Z, and XZ for d = 4. Corre- 
spondingly, there are no other bases which are mutually 
unbiased with respect to these three [13]. 



VI. GENERAL DIMENSION 

Tarry was the first to prove that no two OLSs of or- 
der six exist [20| . i.e., the net for d = 6 has only three 
rows. The operators X m Z n commute for numbers 777 and 
77 from the first cell of these rows and the corresponding 
MUBs are the eigenbases of X, Z, and XZ. Similarly to 
the case of d = 4 no other MUB with respect to these 
three exists [19J . Of course, the question whether differ- 
ent three MUBs can be augmented with additional MUBs 
remains open. 



A. MacNeish's bound 

More generally, the lower bound on the number of 
OLSs was given by MacNeish 21]. If two squares of order 
a are orthogonal, A _L B, and two squares of order b are 



orthogonal, C _L D, then the squares obtained by a direct 
product, of order ab, are also orthogonal, AxC A. BxD. 
This implies that the number of OLSs, C, of the order 
d = p[ x ..-Pr" , with pi being prime factors of d, is at least 
C > mini(p[' — 1), where p^ — 1 is the number of OLSs 
of order p\\ A parallel result holds for MUBs 0,[l9|]. If 
| a) and \b) are the states of two MUBs in dimension d\, 
and |c) and \d) are the states of MUBs in dimension c?2, 
then the tensor product bases \a) ® |c) and \b) <g> \d) form 
MUBs in dimension d\<fo,. Thus, for d = p^ 1 ...p£ n there 
are at least min^pp + 1) MUBs. 



B. Latin operator basis 

In general, we know more about the number of OLSs 
than about the number of MUBs [la]. We use this 
knowledge to derive conditions which restrict the form 
of MUBs. Consider the operators 



d d-l 



B TI 



= i+EE« 



(ii) 



m=0 £=1 



where n m = 0, ..., d - 1 and Sf n = J2j=o rfd \J)m(j\ nave 
complete set of MUBs as eigenbases, m = 0, ...,d. We 
show that existence of such a set and orthogonality of 
d 2 operators -B no ...n d implies completeness of the set of 
OLSs. The trace scalar product Tr(S^ nd B n i ..„/) is 
given by d 2 (k — 1), where k denotes the sum of Kro- 
necker deltas, k = 5„ an > a + ... + 5 nd n' d - Operators B„ „, nd 
and B n i n / are orthogonal if and only if k = 1, i.e. 
77 m = n' m for exactly one 777. This condition applied to 
d 2 orthogonal operators, defines a complete set of orthog- 
onal squares. To see this, take d 2 orthogonal operators 



B 



n (b)...n d (b) 



with b = 1, ..., d 2 and consider d+ 1 squares 
defined by their indices n m (b) for a fixed 777. If the squares 
were not orthogonal, one could find at least two identical 
pairs, (n m (b),n m >(b)) = (n m (b'),n m >(b')), implying that 
operators (fTTj) are not orthogonal (k > 1). Therefore, 
e.g., for d = 6, there is no complete set of MUBs for 
which operators B no nd are orthogonal because there is 
no complete set of OLSs in this case. 



C. Orthogonal functions 

The second condition is obtained by noting that a net 
defines "orthogonal" functions, F a (m,n)^ which give the 
column of the ath row where the pair 777 77 is entered. 
The orthogonality means that for the pairs m 77 for which 
the function F a (m,n) has a fixed value, the function 
F a ' (777, 77) acquires all its values. We show that if d 2 uni- 
taries, U mn , shift (up to a phase) the states of different 
bases in accordance with the net 



Umn\j)a OC \j + F a (m, n)) a , 



(12) 



then these bases are MUBs. For the proof, note that 

Sj'=o U(*K')a'| 2 = 1- From orthogonality of the func- 
tions, this sum can be written as Es \a(j + F a {m,n)\j' + 
F a t(m, n)) a '\ 2 — 1, where S is the set of pairs mn 
for which F a (m,n) has a fixed value. By (|12p. the 
last is SsloO'l^rin^nnliOa'l 2 ) which due to unitarity, 
Umrftmn = 1, is the sum of d identical terms | a (j|j')a'| 2 - 
Therefore, | a (j|j')a'| 2 — V^- Further, given (i 2 unitaries 
with property (|12| , one recovers the table in the following 
experiment: Prepare |0) o , act on it with U mn , measure 
in the same basis, and write the pair m n in the ath row 
and the column corresponding to the result. Thus, in di- 
mension six, there cannot be 36 unitaries satisfying (|12p . 
with the orthogonal functions, for more than three bases, 
because otherwise one could construct more than three 
orthogonal squares of order six, which is impossible. 



VII. HIDDEN- VARIABLE SIMULATION OF 
MUBS 

The net designs can be used to construct hidden- 
variable models which simulate results of complemen- 
tary measurements on certain states. Recently, Spekkens 
showed that only four "ontic states" (hidden variables) 
are sufficient to simulate complementary measurements 
of a qubit prepared in a state of a MUB [22|]. In his 
model, quantum states of MUBs correspond to the "epis- 
temic states" satisfying the knowledge balance principle: 
The amount of knowledge one possesses about the on- 
tic state is equal to the amount of knowledge one lacks 
[22| . This principle lies behind the net design. Left table 
of ([2]) corresponds to the original Spekkens' model: The 
numbers enumerate ontic states, cells correspond to the 
epistemic states and rows to the complementary mea- 
surements. All other tables generalize the model. To 
identify the ontic state one needs two dits of information 
(there are d 2 ontic states) , whereas the epistemic state is 
defined by a single dit, leaving the other one unknown. 
The quantum states described by these models require (a 
classical mixture of) only two dits to model d outcomes 
of d + 1 quantum complementary measurements. 

Our approach allows us to ask the question how many 
epistemic states satisfying the knowledge balance prin- 
ciple, i.e., having d underlying ontic states, correspond 
to quantum states. For example, in the case of a two- 
level system there are four ontic states, and six pos- 
sible epistemic states [see the net design of (J2J]. All 
six correspond to quantum eigenstates of complementary 
observables. In general, any epistemic state is repre- 
sented by a cell of d numbers [i\ 12 ■■■ id]- Since each 
number takes on one of d 2 values, the numbers can- 
not repeat and their order is not important, there are 
£ d = E^iE^+r-'E^tt^+i possible epistemic 
states, with D = d 2 — d + 1. 

For d being a power of a prime the quantum states cor- 
responding to the cells of the net design are basis vectors 



of a complete set of MUBs. They can be used to uniquely 
decompose arbitrary Hermitian operator 



d d-\ 



6 = -Tr(6)i+^^pf ) | J ) m ( J |, 



(13) 



m=0 j=0 



Jrn) 



where p'" 1 ' 1 — m (j\0\j) m and \j) m is the jth state of 
the mth MUB. For the proof, note that the complete 
set of MUBs can be used to define the operator ba- 
sis in the Hilbert-Schmidt space S^ = E 7 =o Vd \j)m(j\- 
There are d 2 such operators, because m — 0,...,d, the 
power £ = 0, ..., d — 1 and all d operators S^ are equal 
to the identity operator. Since they are normalized as 
Tr[(5|j) Sm'] ~ d8 mm i8^ any operator has a unique ex- 
pansion = ±[Tr(6)l + £l=o £ e t! Tr(6(5«jt)£«J 

Writing S^ in terms of projectors on MUBs one finds 
Eq. C3J. 

If O is a quantum state, Tr(O) = 1 and Pj 's are prob- 
abilities to observe outcomes related to suitable states of 
MUBs. We consider general epistemic states, not nec- 
essarily those corresponding to the cells of the net de- 
sign. Such epistemic states have "partial overlap" with 
the cells, defined as the number of common ontic states 
divided by d. For example, the epistemic state [00 01 20] 
has an overlap of | and 4 with the first and third epis- 
temic state of the first row of table (|4]) , respectively. To 
construct operator O associated with a general epistemic 
state, we take these overlaps to define the probabilities 
p\ ■ Since we would like to see how many epistemic 
states correspond to quantum states we take operators 
6 with a unit trace. If Tr(d 2 ) = 1 and Tr(0 3 ) f 1, 
the operator O cannot represent a quantum state, be- 
cause the first condition excludes mixed states, and both 
of them exclude pure states [23] . We find that for d = 3 
only the epistemic states of the net design correspond to 
the quantum states. There are Q3 = 12 such states, 
out of £3 = 84 different epistemic states. The ratio 
of IZd = Qd/£d rapidly decreases with d: we checked 
TZ 3 = 1/7, K A = 8/455 and TZ 5 = 1/1771. Thus, 
most of the epistemic states, constructed according to 
the "knowledge balance principle," do not represent a 
quantum-physical state. 



VIII. CONCLUSIONS 

In conclusion, we showed a one-to-one relation between 
OLSs and MUBs, if d is a power of a prime. For general 
dimensions, we derive conditions which limit the struc- 
ture of the complete set of MUBs and we presented par- 
allelism between the MacNeish's bound on the minimal 
number of OLSs and the minimal number of MUBs. In- 
terestingly, the MacNeish's bound is known not to be 
tight. There are at least five OLSs of order 35, where 
the MacNeish's bound is four 24]. Therefore, further in- 
sight into the relations between MUBs and OLSs would 



be gained from studies of MUBs for d — 35. 

Finally, using the squares, we constructed hidden- 
variable models that efficiently simulate measurements of 
MUBs. However, the majority of states in these models 
do not have quantum-physical counterparts. 
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